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Hybrid graphene-superconductor devices have 
attracted much attention since the early days 
of graphene research[l— 18]. So far, these stud¬ 
ies have been limited to the case of diffusive 
transport through graphene with poorly defined 
and modest quality graphene-superconductor in¬ 
terfaces, usually combined with small critical 
magnetic fields of the superconducting electrodes. 
Here we report graphene based Josephson junc¬ 
tions with one-dimensional edge contacts [19] of 
Molybdenum Rhenium. The contacts exhibit 
a well defined, transparent interface to the 
graphene, have a critical magnetic field of 8 Tesla 
at 4 Kelvin and the graphene has a high qual¬ 
ity due to its encapsulation in hexagonal boron 
nitride[19, 20]. This allows us to study and exploit 
graphene Josephson junctions in a new regime, 
characterized by ballistic transport. We find that 
the critical current oscillates with the carrier den¬ 
sity due to phase coherent interference of the 
electrons and holes that carry the supercurrent 
caused by the formation of a Fabry-Perot cavity. 
Furthermore, relatively large supercurrents are 
observed over unprecedented long distances of up 
to 1.5 n m. Finally, in the quantum Hall regime 
we observe broken symmetry states while the 
contacts remain superconducting. These achieve¬ 
ments open up new avenues to exploit the Dirac 
nature of graphene in interaction with the super¬ 
conducting state. 

In particular, the chiral nature of the charge carri¬ 
ers in graphene is predicted to give rise to specular 
Andreev reflection [21], and the conventional quantum 
Hall effect can be markedly different due to the interac¬ 
tion between edge states and the superconductor [22, 23]. 
Such systems also provide a unique way to probe valley- 
polarized edge states [24], topological confinement in bi¬ 
layer graphene [25], the interplay between superconductiv¬ 
ity and quantum confinement or ballistic two-dimensional 
Josephson junctions and their response to phase coherent 
interference effects. 

There are two important prerequisites that must be 
satisfied in order to observe any of these phenomena exper¬ 
imentally. First, the graphene-superconductor interface 


should be transparent and well defined. Secondly, the 
graphene must be of high electronic quality. In addition, 
for some of the above effects, a superconductor with a 
large upper critical field, H c2 , is required. While signifi¬ 
cant technological progress has been made in improving 
the quality of graphene by either suspending graphene [26] 
or encapsulating it in hexagonal boron nitride (hBN) 
[19, 20], the main challenge has been to combine such 
low-scattering graphene with a (large H c2 ) superconduc¬ 
tor. All reports on graphene-superconductor devices to 
date involved superconducting contacts deposited directly 
on the graphene surface, and diffusive transport through 
the device. In addition to the modest electronic quality, 
the use of top contacts leaves ambiguity in where exactly 
Andreev reflection takes place and under what spectral 
conditions. I.e. it is not clear how far electrons travel 
beneath the contact before entering the superconductor. 

To realise high quality graphene-superconductor junc¬ 
tions we encapsulate graphene between two hBN crystals 
using the van der Waals pick-up method [19] . This method 
ensures that the graphene is never in contact with any 
polymer during the stacking and thereafter. Electrical 
contact is made by metal deposition onto areas where the 
stack has been etched through. Unlike earlier work [19], 
where metal deposition is done in a separate lithography 
step, we start by etching only the region to be contacted, 
followed immediately by metal deposition. This has the 
following advantages: (i) our contacts are self-aligned, 
thereby minimizing redundant metal overlap above the 
graphene and reducing the screening of electric and mag¬ 
netic fields and (ii) combining the etching and deposition 
in one step minimizes resist residues at the contact inter¬ 
face, which is necessary for transparent contacts. Instead 
of a normal metal, we sputter an alloy superconductor 
MoRe, which is attractive in several respects. First, MoRe 
is a type-II superconductor with a critical temperature 
T C «8K and an upper critical field H c 2 » 8 T (at 4.2 K), 
which should easily allow for the observation of quantum 
Hall states while the MoRe remains predominantly su¬ 
perconducting. Secondly, it has been shown that MoRe 
makes good electrical contact to carbon-based materials 
such as carbon nanotubes [27]. Considering the fact that 
edge-contact resistance can vary by an order of magnitude 
depending on the choice of metal [19], it is critical to select 
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Figure 1. High-quality hBN-Graphene-hBN devices, a. 

An optical image of device A. A graphene/hBN sandwich 
(blue) is contacted on both sides from the edge with MoRe 
contacts (yellow). The contacts are split further in two, which 
allows a (quasi-) four probe measurement with minimal series 
lead resistance, b-c. A schematic cross-section of the device, 
d. The measured resistance, R , as a function of gate voltage, 
Vgate , at room temperature and at 4.2K. The carrier density, 
n, is extracted from Shubnikov-de-Haas oscillations, e. Differ¬ 
entiated conductance, dG/dV ga te , as a function of gate voltage 
and magnetic field, taken at 40 mK. f The conductance, G, 
as a function of gate voltage at B = 12 T and T = 40 mK, 
showing the symmetry broken states. 


a superconductor which makes good electrical contact to 
graphene. This is particularly important in the context 
of superconductor (S) graphene (G) JJs, where the trans¬ 
parency of the S-G interface directly affects the Andreev 
reflection. Furthermore, unlike surface contacts, such 
one-dimensional edge contacts ensure that the Andreev 
reflection occurs at a well-defined location, at the edge of 
the graphene, where it contacts with the 3-dimensional 
bulk superconductor. After the deposition of the super¬ 
conducting electrodes, we etch the stack into the desired 
geometry. 

An optical image and a cross-sectional schematic of 
device A are shown in Fig. la-c. The graphene is etched 
to a L = 1.5 fi m long and W = 2.0 pm wide rectangle, with 
MoRe edge contacts on either side. All measurements 
described here are performed in a (dc) four point geome¬ 
try, as shown in Fig. la. The MoRe leads are arranged 
such that the lead series resistance is minimized and the 


measured resistance is effectively the two-probe graphene 
resistance, irrespective of whether the MoRe is normal 
or superconducting. This is important since disordered 
superconductors such as MoRe have a large normal-state 
resistivity, potentially confusing the interpretation of the 
measurements when the electrodes turn normal (see SI). 

Fig. Id shows the measured resistance, i7, versus back 
gate voltage, V ga t e , at room temperature and 4.2 K. A 
clear electron-hole asymmetry is visible with the resistance 
in the hole doped (p) regime being somewhat larger than 
that in the electron doped (n) regime. We attribute this 
to contact-induced n-type doping, which leads to the 
formation of pn junctions close to the contacts when the 
bulk of the graphene is p doped. Such n-type doping 
effects from normal edge contacts have also recently been 
reported [28]. Fig. le shows the Landau fan diagram 
recorded up to B = 12 T. The high electronic quality of 
the graphene is evident from the emergence of broken 
symmetry states above B - 5 T, which are well developed 
at B = 12 T (Fig. If). To our knowledge, this is the first 
observation of broken symmetry states in graphene with 
superconducting contacts. The plateaus on the electron 
side are better developed than those on the hole side, 
presumably a consequence of doping near the contacts. 

At zero magnetic field we observe a gate-tunable su¬ 
percurrent through the device. In Fig. 2a we plot the 
diffential resistance, dV/dl , as a function of gate voltage, 
Vgate , and the current bias, J^ c . Evidently, the critical 
current, J c , vanishes at the charge neutrality point, but 
reaches values in excess of 100 nA at V ga te = 30 V. The 
individual J^ c - V curves are hysteretic, as is evident from 



Figure 2. Long distance Josephson current in edge- 
contacted graphene, a The differential resistance, dV/dl , 
is plotted as a function of applied dc current bias, Ida and 
gate voltage, V ga ta at 60 mK. b Critical current density, J, 
plotted as a function of device length, L. Squares are the side 
contacted MoRe graphene devices A - E reported here. Black 
(red) squares correspond to a temperature of 50 mK (700 mK). 
More details about the temperature dependence can be found 
in the SI. Circles are data points taken from the literature 
[1—18]. Colors indicate different superconductors used: Black 
circles refer to Al, green circles to Nb/NbN/NbTiN, blue ReW, 
red Pb/Pbln and yellow Pt/Ta. 
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Figure 3. Fabry Perot resonances in a Josephson junction, a. The differential resistance, dV/dl , is plotted as a function 
of the applied dc current bias, Idc, and gate voltage, V ga te, at T = 550 mK. At 60 mK, statistical fluctuations in I c make 
the effect much less visible. Inset: A schematic of a cavity formed between pn-junctions due to doping near the contacts. 
Interference occurs due to reflections at the pn-junctions. b. The normal state conductance, Gn, and the critical current, / c , 
plotted as a function of gate voltage, V ga te • c. (upper panel) The conductance measured as a function of the magnetic field and 
gate voltage with Id c = 100 nA. The dispersion of the Fabry-Perot interferences follows a 4 th -order polynomial, see equation 1, 
plotted in yellow, (lower panel) The simulated conductance for a cavity size of L = 1.3 pm and W = 2 pm. d.(upper panel) The 
conductance SG in the npn regime as function of absolute wavenumber |fcp| in the central part of the device (|&f| is determined 
from the carrier density). SGn is obtained after subtraction of a slowly varying background (see SI for details). In the lower 
panel we plot SGn in the nn n regime for the same wavenumber range. Here we attribute the fluctuations to UCF. 


the asymmetry about /^ c = 0 (we discuss the possible 
origins of this hysteresis in the SI). On the hole side I c is 
considerably smaller, consistent with the formation of the 
conjectured npn junctions. In Fig. 2b we plot the critical 
current density per unit length, J, versus the JJ length, 
L, with data obtained from previous reports of graphene 
JJs (circles) along with the present MoRe edge-contacted 
devices (squares). The black squares show the critical 
current density at 50 mK, whereas the red squares are 
taken at 700 mK. We point out that the critical current 
density depends on the temperature and the graphene car¬ 
rier density, which vary from study to study. Despite this, 
it is clear that our MoRe edge-contacted devices stand 
out in relative magnitude compared to the previous data. 
We find large supercurrent densities (up to ~ 200 nA/pm) 
over significantly longer distances (~ 1.5 pm). The obser¬ 
vation of large supercurrents over an unprecedented long 
distance of 1.5 pm indicates a high quality of both the 
graphene itself and of the ID graphene-superconductor 
interfaces. 

In addition, we find unambiguous signatures of ballistic 
Josephson transport in this 2D geometry. As shown in 
Fig. 3a, we observe for the first time clear oscillations in 
the critical current and the retrapping current when we 
vary the gate voltage, indicative of Fabry-Perot (FP) in¬ 
terferences in the supercurrent through the junction. The 
transmission probability of electrons and holes that carry 
the supercurrent is the result of interference of trajecto¬ 
ries that travel ballistically from one contact to the other 
with multiple reflections close to or at the edges of the 
graphene flake. As the gate voltage is varied, the Fermi 
wavelength changes, constructive and destructive inter¬ 
ference alternate, leading to modulations in the critical 
current. One may expect the graphene-superconductor 
interfaces to form the walls of the cavity. However, we 


observe the I c oscillations only on the hole-doped side 
and not on the electron-doped side (see Fig. 3d and the 
discussion below). This suggests that in the presence of 
n-doped regions near the MoRe-graphene interface the 
relevant cavity is instead formed by pn junctions near 
the contacts (see inset Fig. 3a). This gives rise to a 
reduced cavity length L c . This length can be directly 
inferred from the period of the oscillations, extracted via 
a Fourier analysis (see SI) of these oscillations over many 
periods. A cavity length of L c = 1.3 pm is found, which 
is smaller than the etched device length (L = 1.5 pm). A 
similar difference between device size and inferred cavity 
length was seen in device D (see SI). This difference may 
arise from screening of the back gate near the contacts in 
combination with the presence of the n-doped regions at 
the MoRe-graphene interfaces in both devices. 

The interpretation of the oscillations in I c in terms of 
FP interference, is further supported by comparing them 
with the oscillations in the normal state conductance, Gat, 
measured at currents just above I c . The oscillations of 
I c with gate voltage clearly match the oscillations in Gat, 
(see Fig. 3b), as expected for Josephson junctions. In 
the case of normal state transport, we can apply a weak 
magnetic field perpendicular to the graphene, to apply 
a Lorentz force to the trajectories of electrons and holes. 
This is expected to give a characteristic shift of the FP res¬ 
onances due to the accumulation of extra field-dependent 
phases. Indeed in the measurements shown in Fig. 3c, we 
find that as B increases the main resonance features shift 
to higher density, following a characteristic dispersion. To 
enhance the visibility, we plot the quantity G su b, which 
was obtained after subtracting a gate dependent (but field 
independent) modulation of the background conductance. 
We compare the data with the results of numerical simu¬ 
lations of the device conductance (see methods for further 
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details) in the ballistic regime and with npn junctions for 
the exact geometry of the measured device (Fig. 3c lower 
panel). Simulation and experiment show an almost iden¬ 
tical dispersion of the FP resonances with magnetic field. 
It is also possible to obtain a semiclassical expression for 
the resonance condition (see SI) by considering all the 
phases accumulated in the n-region of the npn junction: 

L c 1 1 l L\eB 

A F(V gat e) nm+ 2 + 6n m \ h 

with n m a specific integer mode, A F^Vgate) the Fermi 
wavelength which is tuned by the backgate (V ga te ~ 1/A^), 
L c the cavity size, e the electron charge and h Planck’s 
constant. The yellow curves in Fig. 3d are calculated 
using Eq. 1 for modes n m = -121, -120 and show an 
excellent agreement with the measured and simulated 
results. This provides strong evidence that the observed 
oscillations, both in I c and Gw, arise from Fabry-Perot 
interference, which implies phase-coherent ballistic trans¬ 
port. While such oscillations due to FP interference have 
been reported before in a variety of systems including 
high-quality graphene with normal contacts [2, 30], here 
we provide evidence for phase coherent FP interference 
in the supercurrent, which has not been observed before 
in any 2D geometry. 

In order to better understand the microscopic details of 
our device, we compare the conductance in the npn regime 
with that in the nn'n regime (Fig. 3d). Whereas in the 
npn regime (upper panel), we observe periodic oscillations 
as a function of absolute wave number, |fcp|, we observe 
universal conductance fluctuations (UCF) in the nn'n case 
(lower panel). We attribute these fluctuations to diffuse 
boundary scattering at or close to the graphene-MoRe 
interface. This diffuse scattering should also be present 
on the hole-side but does not dominate the transport due 
to the presence of the pn junctions. Using the ballistic 
limit, L much larger than the mean free path, where all 
resistance is from the contact interface, we can estimate 
a lower bound on the contact transparency, T via G = 
From the conductance in the nn'n regime 
(see SI) we find a contact transparency of T > 0.2. In 
the npn case, the conductance is dominated by the pn 
barriers. In this case, we can estimate the sharpness, d, of 

2 fT 

the p to n transition regions via G npn = -y-FF. We 
find a sharpness of d ~ 70 nm, which is a plausible value 
considering the device dimensions. 

Since the DC Josephson effect is observed in these 
graphene devices over micron scale distances, we can also 
explore the magnetic field dependence of the critical cur¬ 
rent for unusual geometries. Earlier reports concerned 
graphene Josephson junctions with lengths much shorter 
than their width. In this case, the magnetic field depen¬ 
dence of I c is expected to follow the standard Fraunhofer 
diffraction pattern observed in tunnel junctions [31]. In 
the present devices, in contrast, the aspect ratio is close 
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Figure 4. Anomalous Fraunhofer diffraction pattern, a. 

The differential resistance ( dV/dl ) is plotted as a function 
of applied current bias, Ida an d magnetic held, £>, at a gate 
voltage of 30 V. We observe a a separation between minima 
that clearly exceeds the flux quantum, h/ 2e. 


to 1, which has two consequences. First, unlike in tunnel 
junctions, the phase difference across the junction must be 
integrated along both interfaces. Furthermore, contribu¬ 
tions involving reflections off the side of the junction must 
be included, especially when transport is ballistic [32-34]. 
The main prediction in this case is that the periodicity 
of I c with magnetic flux becomes larger than a single 
flux quantum, 4>o = h/ 2e. Despite significant differences 
across the patterns measured on the various devices, we 
consistently find a period larger than 4>o, as seen in Fig. 4 
for device A (and in the SI for devices B and G). In con¬ 
trast, earlier reports on graphene Josephson junctions all 
show flux periods smaller than 4>o [1, 2, 5, 7, 9, 12, 13, 16] 
before corrections to account for the London penetration 
depth. 

The Fabry Perot oscillations in the critical current and 
the anomalous Fraunhofer diffraction patterns are mutu¬ 
ally consistent and provide strong evidence of ballistic 
effects in superconducting transport through graphene. 
We believe that this is the first unambiguous demonstra¬ 
tion of a ballistic JJ in graphene. 


Methods 


DC transport measurements 

Low temperature dc measurements are done in a Leiden 
Cryogenics MCK-50 3 He/ 4 He dilution fridge. The setup 
can reach a base temperature of 40 mK with an electron 
temperature of about 70 mK. DC Currents and voltages 
are applied and probed with a home-built measurement 
setup. Furthermore the setup is equipped with a super¬ 
conducting magnet coil that can sustain fields of up to 
12 T. 
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Tight-binding simulation 

The FP oscillations in the npn junction are simulated 
by a tight-binding calculation using the Kwant software 
package [35]. The source code has been provided along 
with this submission as an ancillary file. A 1.5 /imx2.0 pm 
hexagonal lattice is discretized with a lattice constant 
of a = 2 nm, with metallic leads on the 2.0 pm wide 
sides. The contact induced doping near both leads is 
modeled by a 100 nm region with a fixed chemical po¬ 
tential. The width of the transition region from the n 
to the central p region is set to 50 nm and modelled by 
tanh [(x - xo)/25 nm]. A finite contact resistance is im¬ 
posed by reducing the transparency between the central 
strip and the leads to 60%. Finally we calculate the trans¬ 
mission as a function of the Fermi wavenumber kp(p v ) 
and magnetic field B, resulting in the dispersion given in 
Fig. 3d. 
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SUPPLEMENTARY INFORMATION 

DEVICE FABRICATION 

Graphene flakes are mechanically exfoliated from natu¬ 
ral graphite (NGS naturgraphit GmbH) with blue Nitto 
tape (SWT 20+ Nitto Dekko Corp.) and put on Si sub¬ 
strates with a 90 nm thick thermally grown Si 02 top 
layer (IDB technologies Ltd.). Single layer graphene is 
located under an optical microscope (Olympus BX51). 
The bottom hexagonal Boron Nitride (hBN) flakes are 
made in the same way as the graphene flakes, using sub¬ 
strates with 285 nm thermally Si 02 (NOVA Electronic 
Materials, LLC). The top hBN flakes are prepared on 
polymer coated transparent substrates. The transpar¬ 
ent substrates consist of a microscope glass slide with a 
~ 1 x 1 cm 2 PDMS (Polydimethylsiloxane) film(GelPak, 
WF-40-X4-A). On top of that we spin (1500 RPM for 55 
sec) 9 wt% MMA (Methylmethacrylate) in Ethyl lactate 
and bake it for 20 min at 120 °C. The MMA film acts as a 
transfer agent and the PDMS film promotes a cushioning 
effect during the pickup and transfer. The glass slide with 
hBN flakes is mounted in a micromanipulator. The setup 
is similar to that of ref. [ 1 ] but with an additional heating 
stage implemented. The graphene on Si/Si 02 substrate 
is heated to about 80 °C and brought in contact with 
the top hBN. The MMA becomes soft and the graphene 
is picked up through adhesion to the hBN. Subsequent 
transfer is done at a higher temperature of 100 - 150 °C. 
The stack is patterned using standard E-beam lithogra¬ 
phy with a thick PMMA mask and is shaped by reactive 
ion etching in a 20 ml min -1 O 2 plasma during 1 min 
(Leybold-Heraeus, 60 W, 50 /ibar) and subsequently in 
a 40/4 ml min -1 CHF 3 /O 2 plasma during 1 min (60 W, 
50 /ibar). Next 80 nm MoRe is deposited by a sputter¬ 
ing process using a DC plasma with a power of 100 W 
in Ar. Lift-off is completed in hot acetone (T = 54 °C) 
for about 3-4 hours. A second lithography step defines 
the intended graphene geometry in a PMMA / hydro¬ 
gen silsesquioxane (HSQ) mask, followed by an etch step 
using the same plasma conditions as in the first lithog¬ 
raphy step. The fabricated devices are not exposed to 
any (Ar/H 2 ) annealing processes that are common for 
graphene devices. 


DETAILS ON THE npn-JUNCTION 
SEMICLASSICS 

In Fig. 3c of the main text we have highlighted the 
dispersion of dominant Fabry Perot resonances with the 
external magnetic field. Here we derive these resonances 
from the semiclassical phase accumulated in the p-region 
of the npn-junction. If the junction is along the re¬ 
direction the total Landauer conductance of the junction 
can be expressed as a sum of transmission amplitudes 


over transverse momenta k y . The leading contribution of 
the oscillating part of the conductance is given by [ 2 ] 

o ^ 2 

Goec = — E |T + | 2 |T_| 2 | J R + || J R_|cos0e- 2L ^. (2) 

h b, 

Ky 


Here L c is the length of the p-doped region (the cavity), £ 
the mean-free path in the p-region. Choosing the center of 
the junction at x - 0 and using the Landau gauge where 
= Bxy, the transmission and reflection amplitudes at the 
junctions located at x - ±L /2 are 


T± = exp 


7 r hv? ( eBL \ 2 

r^ir) 


and |i?±| = + - |T±| 2 . Finally 


( 3 ) 


< 4 > 

is the semiclassical phase difference across the cavity for 
neighboring trajectories. For simplicity we have assumed 
that the junctions at x ± L /2 are sharp. We are inter¬ 
ested in resonances at finite fields where the transmission 
is dominated by small transverse momenta (see below) 
which acquire an additional nontrivial i r phase, up on 
reflection from the pn-junction. The total accumulated 
phase across the cavity and back is 6 - 2 #wkb + tt- At zero 
magnetic field k y - 0 is equivalent to normal incidence 
yielding a perfect transmission (Eq. (3)) and a vanishing 
contribution to the interference. In this regime the dom¬ 
inant contribution to the Fabry-Perot resonances stems 
from finite momenta. The bending of the semiclassical tra¬ 
jectories due to the cyclotron motion at finite field causes 
a switch in the dominant contribution of the conductance 
to small momenta. The transition is characterized by a 
phase shift in the conductance oscillations at [ 2 ] 
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For fields B larger than the characteristic field the 
largest contribution to the total conductance Eq. (2) is 
from k y - 0 (|T + | 2 |T_| 2 |R+||R_| is maximized by k y - 0). 
The switching field dependents on the degree of collima- 
tion at the np/pn-junctions which is determined by the 
electrical field E = V/d, where V = (/i n -/i p )/e across the 
junction. An estimate for the junction thickness d- 70 nm 
is given in the main text. Unfortunately it is not possible 
for us to get a reliable estimate of the doping level in the n- 
doped region near the contacts. We observe pnp-junction 
behavior up to very large negative back gate voltages. 
This suggest that screening from the contacts plays an 
important role in the local density level and profile render¬ 
ing even a crude estimate on simple grounds impossible. 
From our experimental measurements we estimate the 
switching field to be below 20 mT. This implies that the 










local electric field should be on the order of 1V/ fi m or 
smaller. The field range studied in Fig. 3c is well above 
this switching field allowing us to focus on the regime 
where the k y - 0 contribution to the conductance oscilla¬ 
tions dominates. This contribution is resonant when the 
total acquired phase for k y = 0 is an integer multiple of 
27r 

26 > WKB|/c y =0 + 7T = 27m. (6) 

For small fields where the cyclotron radius is larger than 
the junction length L we can expand #wkb to lowest 
order in B 

( 2r2o2\ 

'--wst)- (7) 

Using k F ~ 270/im -1 and B = 100 mT we obtain a cy¬ 
clotron radius r c = hk F /eB = 1.8 /im. This confirms that 


our experiments are still in this regime. Finally from the 
resonance condition equation (6) we find 

k F L = n7r + - + — (L 2 -B ) 2 . ( 8 ) 

2 6n{ h J K J 


Equal contributions 
t L.M.K.Vandersypen@tudelft.nl 

[1] Andres Castellanos-Gomez, Michele Buscema, Rianda 
Molenaar, Vibhor Singh, Laurens Janssen, Herre S. J. 
van der Zant, and Gary A. Steele. Deterministic transfer of 
two-dimensional materials by all-dry viscoelastic stamping. 
2D Materials , 1(1):011002, 2014. 

[2] Andrea F. Young and Philip Kim. Quantum interference 
and Klein tunnelling in graphene heterojunctions. Nat. 
Phys ., 5:222-226, 2009. 




9 



/* ( n A) 


b 2A/ne 



140 

120 

100 

_ 80 
< 

J 60 

40 
20 
0 

0 200 400 600 800 

T (mK) 



Figure SI. Current-voltage characteristics and temperature dependence a. A current-voltage curve (in device A) 
at 50 mK and large n-doping showing a characteristic hysteresis. Such hysteresis is often attributed to the existence of an 
underdamped Josephson junction. An alternative explanation for such hysteresis lies in non-equilibrium or heating effects in the 
JJ. In this scenario increasing the current results in an increased power fed through the JJ, Joule heating, and an increased 
electron temperature compared to the temperature of the reservoirs. Both these processes could be at play in our graphene JJs, 
and one would require further detailed studies to determine the origin of the hysteresis, b. Bias spectroscopy (device E) shows 
no discernable features associated with multiple Andreev reflections (MAR). Top axis shows the expected positions of these 
features. There can in principle be several reasons for the absence of discernable MAR features in the differential resistance, 
such as extremely low (or high) transparency or the presence of residual inelastic scattering. To resolve this question we need 
more data from a range of devices, a more accurate determination of the contact transparency, and a more mature theoretical 
model for ballistic graphene-superconductor devices, c. Temperature dependence of I c and I r t for device A. We can compare 
this data with the behavior predicted by either the Eilenberger equations (for clean SNS junctions) or the Usadel equations (for 
diffusive junctions). Both models predict a smooth variation of I c with temperature, so the pronounced kink around 300 mK 
requires a separate explanation. 



Figure S2. Fabry-Perot interference in device D . a. False colored scanning electron microscope image of device D with 
MoRe contacts. The measurements in panel b are recorded between the contacts indicated in the figure. The distance between 
the contacts is 1 /im and the contact width is 360 nm. b. The conductance G between the contacts indicated in panel a plotted 
as function of magnetic field B and gate voltage. Here we observe a very similar parabolic dispersion as in Fig. 3c in the main 
text. The yellow lines are calculated via equation 1 in the main text for L = 0.8 /im and rim = -56, -55. These results imply 
that the transport across the Hall bar device is ballistic. 
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Figure S3. Power spectral density, a. The normal state conductance Gn measured in the npn regime as a function of 
number &f- The dashed line is the slowly varying background, defined here by a 6 th -order polynomial fit. b. The normal state 
conductance Gn in the nn'n regime as a function of wavenumber &f, with the dashed line the slowly varying background. By 
taking the discrete Fourier transform T{5Gn(J^f)} = T(A) of the conductance fluctuations in a and b (after subtracting the 
steep background, see also Fig. 3d, we can obtain the power spectral density PSD( A) = T(A)T*(A). c. The power spectral 
density PSD for device A in the npn and nn'n regime and for device D in the npn regime as a function of length cc/2. Device 
A and D have a designed length of La = 1.5 /im and Ld = 1.0 /im, respectively. From the peak positions in the PSD we extract 
a cavity size of L^ p = 1.3 /im and L^ p =0.8 /im. We attribute the difference of ~ 200 nm to two effects. Firstly, the size of the 
graphene/hBN is most likely slightly smaller than the designed size due to the etching process. Secondly, the cavity is formed 
between the pn junctions, where the n -region is caused by doping from the contact and has a finite length. This would mean the 
n-doped region is about ~ 100 nm on each side. The fact that the pn junctions act as the cavity walls leading to FP interference, 
is further supported by the fact that PSD in the nn'n regime shows no predominant period. 



d 





0 2 4 6 8 10 12 

fl( T) 


Figure S4. Lead resistance switching at 4.2 K. a. An optical image of the die of device D. In this device, transport is 
measured between contacts A and B in a two-terminal configuration. In the normal state, the MoRe leads have a non-negligible 
resistance. This causes artifacts in the measurement when switching between the superconducting and normal state of the 
MoRe. Such artifacts from the lead resistance are minimized in the quasi four-terminal configuration used in device A, see 
Fig. la in the main text. b. Landau fan diagram of Device D (shown also in Fig S2a), showing the conductance as a function 
of magnetic held B and gate voltage V ga te- Above the critical held of MoRe H c 2 ~8 T, the MoRe leads switch to the normal 
state, causing a jump in the two-terminal conductance, c. The conductance G in units of e 2 //i for filling factors v = 2, 6 and 10 
as a function of magnetic held B. d. When plotting resistance instead of conductance, it becomes clear that there is a hxed 
resistance increase when the MoRe leads turn normal. To facilitate comparison, we plot R — where R u = is the ideally 
expected resistance for filling factor v. e. The square resistance R n of the MoRe sheet measured using a Hall bar as a function 
of the magnetic held B. We observe again the upper critical held of MoRe of about H c 2 ~8 T 
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Figure S5. Anomalous Fraunhofer diffraction patterns, a. The differential resistance ( dV/dl ) as a function of applied 
current bias, /d c , and magnetic field, B at a gate voltage of 30 V, for device A. This is the same plot as in the main text, Fig. 4. 
b-c. Similar plots for devices B and C. d. The same measurement as in c, but shown over a larger magnetic field range. The 
lobes persist up to at least 20 mT which corresponds to about 9-10 periods. In general, the dependence of I c on B shows a wide 
range of behaviors, seen through the various Fraunhofer diffraction patterns in the four panels. An important question is what 
kind of variations in our samples contribute to the variation in the Fraunhofer response. For instance, we do not know a priori 
whether the edge contacts have uniform transmission over the full width of the contact. Despite the variations, we consistently 
observe that (1) when deducing the flux periodicity from the first lobe, the flux periodicity is larger than one flux quantum 4>o 
and (2) the minima in I c do not reach zero between the lobes. 















